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We develop a model inspired by the Auger mass composition for Ultra High Energy Cosmic Rays (UHECR). 
The crucial experimental fact for the model is given by the strong proton dominance in energy range (1-3) EeV in 
agreement with the HiRes data. At higher energies the Auger data show progressively heavier mass composition 
reaching the Iron-dominated one at i? ~ 35 EeV. Assuming extragalactic origin of this EeV proton component, 
we argue that it disappears at higher energies due to low maximum-energy of acceleration for the protons -E^ax ~ 
(4 — 10) EeV. For rigidity acceleration mechanism the maximum energy for nuclei is Z times higher, where Z 
is a nucleus charge number, and thus maximum energy for Iron nuclei does not exceed (100 - 300) EeV. This 
energy is close to the maximum energies observed by Auger. The cutoff of the spectrum is provided by photo- 
disintegration of the nuclei strengthened by acceleration cutoff. This model has disappointing consequences for 
future observations in UHECR: (i) Since the energies per nucleon are less than (2 - 5) EeV, the photopion 
processes in extragalactic space are absent together with GZK cutoff and production of cosmogenic neutrinos; (ii) 
Correlation with the nearby sources even at highest energies is absent, because of the nuclei deflection in galactic 
magnetic fields. 



1. Introduction 

There is a dramatic conflict between the re- 
cent observational data of two largest UHECR 
detectors: HiRes [T] and Auger [5]. The HiRes 
data confirm well the signatures of proton propa- 
gation through CMB, GZK cutoff [3] and pair- 
production dip [5 - [H], together with proton- 
dominated mass composition [S] . The Auger data 
strongly favour the nuclei composition progres- 
sively heavier in the energy region (4 - 40) EeV 
and indicate the strong spectrum steepening at 
the highest energies not much consistent with the 
predicted shape of the GZK cutoff. 

We shall discuss first the HiRes data. To study 
the signatures of proton propagation through cos- 
mic microwave background radiation (CMB) we 
assume that generation spectrum of protons is 
proportional to E~'^!> and sources are homoge- 
neously distributed over the universe. Propa- 
gating through CMB protons obtain the spec- 
tral feature in the form of the pair-production 
dip due to p -I- 7cmb p + e'^ -\- e~ scattering 
and the GZK cutoff due to photopion production 
P + 7cMB — > iV + TT. It is convenient to study 



the theoretical energy spectra in the form of the 
modification factor ri{E) given as a ratio of en- 
ergy spectrum Jp{E) calculated with all energy 
losses taken into account, and unmodified spec- 
trum Jp"™, where only adiabatic energy losses 
(red shift) are included. 

r^{E) = J,{E)/j;-^{E). (1) 

Modification factor is less model-dependent quan- 
tity than the spectrum. In particular, it depends 
weakly on generation index 7^, because both 
numerator and denominator in Eq. ([1]) include 
E^'^" . The dip and beginning of the GZK cutoff 
in terms of the modification factor do not depend 
of distances between sources, different modes of 
proton propagation (from rectilinear to diffusion) , 
local ovcrdensity and deficit of the sources etc (see 
[5]). The modification factor is presented in com- 
parison with HiRes spectrum in the left panel of 
Fig. [TJ If to switch off all interactions except the 
adiabatic, one obtains the horizontal line ?/ = 1 
shown in the plot. If to include additionally pair- 
production, one obtains the dip shown in the fig- 
ure by curve rj^f.. The photopion energy losses 
result in the steepening of the spectrum (GZK 
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Figure 1. Left panel: Pair-production dip and GZK cutoff in terms of modification factor in comparison 
with the HiRes observational data. Curves jytot and r]ee show the total spectrum and spectrum with 
adiabatic energy losses and pair-production, respectively. Right panel: E1/2 as numerical characteristic 
of the GZK cutoff in the integral HiRes spectrum (see text). 



cutoff), which starts at energy E ^ 'bO EeV. 
The theoretical dip in Fig. [T] has two flatten- 
ings. The low-energy flattening at i? 1 EeV 
provides transition to galactic cosmic rays, since 
the steep galactic component (oc E~^'^) unavoid- 
ably intersects the flat extragalactic spectrum at 
E <\ EeV. The high-energy flattening explains 
the ankle observed at ~ 1 x 10^^ eV. Compari- 
son with HiRes data shows that both signatures of 
proton interaction with CMB are well confirmed: 
the dip is seen with very good x^i while steepen- 
ing of the spectrum which starts at i? ~ 50 EeV 
agrees within limited statistics with the theoreti- 
cal shape of the GZK cutoff. The nature of this 
steepening as the GZK cutoff is further confirmed 
by the right panel of Fig. [T] valid for the integral 
spectrum. 

In the integral spectrum the GZK cutoff is 
characterized by energy i?i/2, where calculated 
spectrum J(> E) becomes half of the power- 
law extrapolation spectrum KE~"' from low en- 
ergies. As calculations ^ show, this energy is 
Ei/2 = 10^^'^^ eV for a wide range of generation 
indices from 2.1 to 2.8. HiRes collaboration found 
Ei/2 = 10^^-^'^='=° °^ eV in a good agreement with 
the theoretical prediction. In the right panel of 
Fig. [1] we reproduce the HiRes graph [10] from 



which Ei/2 was determined. The plotted value is 
given by ratio of measured flux J(> E) and its 
power-law approximation KE~^ . Extrapolation 
of this ratio to the higher energies is given by 
unity, while intersection of measured ratio with 
horizontal line 1/2 gives £'1/2- 

With some caution one may conclude that 
HiRes has detected the signatures of proton inter- 
action with CMB in the form of pair-production 
dip and the GZK cutoff. The final proof of this 
conclusion must come from direct measurement 
of the mass composition of primaries, and HiRes 
data on fluorescent measurements of elongation 
rate and distribution over Xmax, the atmospheric 
height of the shower maximum, confirm the dom- 
inance of proton composition, indeed. In Fig. [2] 
we plot the data of HiRes [3] on Xmax as a func- 
tion of energy (elongation rate) and RMS(Xi„ax), 
the width of distribution over Xmax- One can see 
that both quantities agree with proton-dominated 
composition. 

Auger data on spectra and mass composition 
are quite different. In contrast to HiRes, Auger 
data [11] show nuclei mass composition starting 
with E ^ A EeV, which becomes progressively 
more heavier with increasing energy (see Fig. [3]). 
The energy spectrum has a sharp steepening at 
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Figure 2. X^ax as function of energy (left panel) and RMS(Xmax), the width of distribution over X,„ax, 
(right panel) from HiRes data |9]. The calculated values for protons and Iron are given according to 
GGSJETII model [E]. 



E (30 - 40) EeV, but energy shape of this 
steepening is quite different from one predicted 
for GZK cutoff. The width of Xmax distribu- 
tion is a very powerful tool for determination of 
mass composition, free of many uncertainties in- 
volved in method of measuring the absolute value 
of Xinax- The narrow width found in Auger ex- 
periment is difficult to falsify, but on the other 
hand the whole picture obtained by HiRes looks 
self-consistent. In particular, the confirmation of 
pair-production dip in four experiments, includ- 
ing Auger |5|8j is a strong experimental argument 
in favour of proton-dominated composition. 

However, in this paper we address the question 
what consequences follow from the mass com- 
position and energy spectrum measured only by 
Auger detector, and what is most natural model 
for explanation of these data. 

2. The model 

We consider here a physical model which di- 
rectly follows from the Auger observations. Some 
additional model ingredients are accepted as as- 
sumptions, which will be justified in the end of 
this section. 

The experimental fact, basic for our model, is 



the observed proton composition in energy range 
(1 - 3) EeV, which is common for both Auger 
and HiRes (see Fig. [31 especially left panel for 
Auger data and Fig. [2] for HiRes data). We 
combine this observation with assumption that 
these protons are extragalactic. The third, and 
the last ingredient of our model, is assumption 
of rigidity-dependent acceleration in the sources, 
formulated in terms of maximum acceleration en- 
ergy as i?max = ^^0, where Eq is an universal 
energy to be determined from the data and Z is 
atomic charge number. These three assumptions 
complete the definition of the model. 

We determine first the maximum acceleration 
energy for protons E^'^^ — Eq . For this we calcu- 
late the extragalactic diffuse proton flux assum- 
ing the power-law generation spectrum Qg(E) oc 
E'''!' with -Emax = Eq and normalizing the cal- 
culated flux by the Auger flux at (1 - 3) EeV. 
For 7g in a range 2.0 - 2.8 we find the maximum 
value of Eq allowed by Auger mass composition 
and spectrum (increasing Eq beyond this limit 
one has a contradiction with mass composition 
and spectrum in Fig. [3]) . 

The results are presented in Fig. U) In cal- 
culations we use a homogeneous distribution of 
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Figure 3. Auger data [TT] on X^ax as function of energy (left panel) and on RMS(Xinax), the width of 
distribution over Xmax, (right panel). The calculated values for protons and Iron are given according to 
GGSJETII model [12]. One can see from the right panel that RMS distribution becomes more narrow 
with increasing energy which implies the heavier composition. 



the sources without cosmological evolution (evo- 
lution parameter m = 0) and with maximal red- 
shift Zmax = 4. As a criterion of contradiction 
we choose an excess of calculated proton flux 
at energy (4 - 5) EeV, where Auger data show 
the dominance of nuclei. The contradiction has 
a different character for different values of jg. 
For the steep spectra 7g w 2.6 — 2.7 the Auger 
spectrum shape and flux are well described by 
^max ^ ^q20 _ iq21 g^^^ contradiction occurs 
only due to mass composition. The extreme case 
is given by 7g = 2.8 and is displayed in the left 
panel of Fig. 01 All curves with E'max > 10 EeV 
are below the data points a.t E > 5 EeV and 
these curves are excluded only by prediction of 
pure proton composition a.t E ^ (4 — 5) EeV, i.e. 
due to contradiction in mass composition with 
Auger in very narrow energy range. For the flat 
generation spectra the contradiction is very pro- 
nounced: the predicted total proton flux exceeds 
the observed one. It can be seen in Fig. |4] for ex- 
treme case 7^ = 2.0 when for i?max — 5 EeV the 
calculated proton flux exceeds the observed one 
at _E w 2 EeV. We conclude that with some re- 
dundancy E™'^^ ~ (4 — 10) EeV for all generation 



indices in the range 2.0 - 2.8. Maximum energy of 
acceleration for Iron is iJ^^^'^ - (100 - 300) EeV. 

In Fig. [5] we plot, as an example, the total 
UHECR spectrum in the 'disappointing model' 
using 7g = 2.3, which might be the case for accel- 
eration by relativistic shocks. The proton spec- 
trum is calculated here in the diffusive model, 
more realistic for energies below 1 EeV. We as- 
sume the turbulent magnetic field with basic 
scales {Bc,lc) =(1 nG, 1 Mpc), the distance be- 
tween sources 40 Mpc and the Kolmogorov diffu- 
sion coefficient (for notation and method of calcu- 
lations see references [13]). The analysis of proton 
maximum energy of acceleration (see left panel of 
Fig- El ) gives Eq = E^'^'^ = 4 EeV, in a rough 
agreement with analysis for homogeneous distri- 
bution of the sources. In fact this set of param- 
eters gives the lowest Eq allowed by Auger data, 
because prediction for maximum energy for Iron 
nuclei is 100 EeV, while in Auger particles with 
almost 200 EeV are observed. However, changing 
the parameters, in particular increasing 7^, it is 
easy to increase Eq to (5 - 6) EeV. 

The inclusion of diffusion provides ffattening 
of proton spectrum below 1 EeV, seen in Fig. [5] 
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Figure 4. Comparison of calculated proton spectra with the combined Auger spectrum for different E™'^^. 
The cases 7^ = 2.8 and 2.0 are shown in left and right panels, respectively. The allowed i^max correspond 
to curves below experimental point aX E ^ A EeV, from which nuclei start dominating. 



as 'diffusive cutoff', because flux J{E) is multi- 
plied to E^ . The 'diffusive cutoff' provides the 
transition from a steep galactic spectrum, most 
probably composed by Iron, to a flat spectrum of 
extragalactic protons. 

The spectrum of nuclei in Fig. [5] is obtained 
by subtraction procedure first suggested in [TJ. 
We subtract from the total Auger spectrum the 
proton spectrum as calculated above, and the re- 
sultant flux is plotted in the right panel of Fig. [5] 
as sum of different nuclei species. In the mixed 
composition model [IS] with arbitrary fractions 
of all primary nuclei at generation it looks easy 
to fit the obtained nuclei spectrum. 

However, the basic feature of the Auger mass 
composition - progressively heavier composition 
with increasing energy - is guaranteed by rigidity- 
dependent maximum energy of acceleration in our 
model: at energy higher than ZEq the nuclei with 
Z disappear, while more heavier (with larger Z) 
survive. Starting from iS™'''' (4 - 10) EcV, the 
higher energies are reachable only for nuclei with 
progressively larger Z. In particular, the max- 
imum observed energy must correspond to Iron 
nuclei, which can reach (100 - 300) EeV. We plan 
to perform the detailed calculations in our next 
paper. 

The predictions of our model are very disap- 



pointing for future detectors. 

The maximum energy, (100 - 300) EeV for 
Iron nuclei, implies only (2 - 5) EeV per nucleon. 
Therefore, practically no photopion processes are 
induced by UHECR and no fluxes of cosmogenic 
neutrinos are produced. 

Correlation with UHECR sources are absent 
because nuclei are deflected in galactic magnetic 
field. The absence of the correlation is strength- 
ened in our model by dependence of the maximum 
energy on Z . 

The signatures of 'disappointing model', which 
can be found by Auger detector, are mass-energy 
relation, seen as elongation rate X^ay^{E), and 
transition from galactic to extragalactic cosmic 
rays below the characteristic energy ~ 1 EeV. 

There are some uncertainties in the calculations 
above. Most important ones are related to esti- 
mates of E^'^^. The model collapses if this energy 
increases to e.g. (50 - 100) EeV. 

Due to uncertainties in Auger data the mass 
composition at (1-3) EeV can include some ad- 
mixture of nuclei, and proton flux is lower. This 
effect decreases E^'^^, and thus is helpful. 

Another effect which slightly decreases E^^^ is 
cosmological evolution of the sources, and it was 
the reason we did not include it in calculations 
shown in Fig. H) 
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Figure 5. Left panel: Comparison of calculated proton spectra with the combined Auger spectrum for 
7g = 2.3 and diffusive proton propagation (see text for details). The cutoff at Ep^^^ — 4 EeV is needed 
not to contradict the data at £^ > 3 EeV. Right panel: Total UHECR spectrum in 'disappointing model' 
in comparison with the combined Auger spectrum. Spectrum of protons is taken from the left panel. The 
spectrum of nuclei is obtained by subtraction procedure as described in text. 



EeV protons detected by Auger can be sec- 
ondary, i.e. produced by photo-dissociation of the 
primary nuclei. In |18| it is demonstrated that 
flux of secondary protons in EeV range is always 
smaller than flux of parent primary nuclei, and in 
[19] considerably smaller than the sum of primary 
and secondary nuclei. 

E^'^^ is calculated for homogeneous distribu- 
tion of the sources, when all modes of propagation 
result according to propagation theorem [20] in 
the same universal spectrum. In case of diffusive 
propagation, for large enough distance between 
sources, E'p'^^ decreases with all results remain- 
ing intact. 

3. Conclusions 

The model we suggested here is aimed to ex- 
planation of observational data of the Auger de- 
tector only. The observational feature crucial for 
the model is the proton composition in energy 
range (1 - 3) EeV. There are two main assump- 
tions of the model: observed protons are extra- 
galactic and they are accelerated by rigidity de- 
pendent mechanism with i?max = ZEq, where Eq 
is the same for all nuclei. The upper limit on 



Eq (maximum acceleration energy for protons) 
is obtained calculating the proton spectrum at 
higher energies using the generation index 7g and 
normalizing flux at (1 - 3) EeV by Auger data. 
Cutoff of the proton spectrum at Eq must pro- 
vide absence of contradiction with Auger flux and 
mass composition at higher energies. The calcula- 
tions are performed for homogeneous distribution 
of sources without assumption of their cosmolog- 
ical evolution, with maximum redshift Zmax = 4 
and for the wide range of 7^ from 2.0 to 2.8. The 
obtained upper limit on Eq is (4 - 10) EeV. The 
maximum predicted energy corresponds to Iron 
and equals to (100 - 300) EeV. The maximum 
energy per nucleon is only (2 - 5) EeV, and photo- 
pion processes on CMB radiation are practically 
absent. Therefore, the GZK cutoff and produc- 
tion of cosmogenic neutrinos are absent, too. The 
cutoff of the spectrum is provided by photodisin- 
tegration and is strengthened by acceleration cut- 
off. 

The rigidity-dependent Smax provides energy- 
dependent mass composition: at energy higher 
than ZEq the nuclei with Z disappear, while more 
heavier (with larger Z) survive. It agrees qualita- 
tively with Auger observations. This feature dis- 
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favours the correlation with the UHECR sources 
at the highest energies: at 100 EeV the Iron nuclei 
dominate and their deflection in galactic magnetic 
fields prevents the correlation with the sources. 

There are two signatures of the model, which 
Auger can observe. The realistic one is transition 
from galactic to extragalactic cosmic rays. This 
transition occurs due to intersection of flat extra- 
galactic proton spectrum with steep spectrum of 
galactic iron nuclei. The flat spectrum of protons 
below 1 EeV appears due to diffusion of protons 
in extragalactic magnetic field [16|17j . The tran- 
sition below 1 EeV combined with nuclei compo- 
sition above (3 - 5) EeV is a strong signature of 
the discussed model. 

The second signature is energy-dependent mass 
composition, heavier with increasing the energy, 
indication to which is already seen in the Auger 
data. 

Is there any alternative explanation of EeV pro- 
tons in the Auger data? 

The EeV protons could have a galactic ori- 
gin, though it contradicts the Standard Model for 
Galactic Cosmic Rays (see e.g. jll]), where max- 
imum acceleration energy, 1 x 10"'^^ eV is reached 
by Iron nuclei. However, one may assume an ad- 
ditional high-energy-proton component of galac- 
tic cosmic rays extended to (2 - 3) EeV. The usual 
argument with anisotropy can be by-passed con- 
sidering a case suggested in [23]. GRB could oc- 
cur in our galaxy 10^ — lO'' yr ago producing high 
energy protons by inner shocks acceleration. The 
protons propagate in the mode of non-stationary 
diffusion, so that most of particles have already 
escaped from the galaxy and now only the tail of 
retarded particles with reduced anisotropy is ob- 
served. In this model the transition from galac- 
tic to extragalactic cosmic rays occurs at ankle, 
which is clearly seen in Auger data. This case 
differs from "standard" ankle model, where tran- 
sition occurs from galactic Iron to extragalactic 
protons. This model needs more detailed investi- 
gation. 

In conclusion, we face at present most serious 
disagreement in observational data of two biggest 
experiments in UHECR. 

HiRes observes signatures of proton propa- 
gation through CMB in the form of the pair- 



production dip and GZK cutoff and these ob- 
servations are well confirmed by direct measure- 
ments of proton-dominated mass composition. 
These observations predict the detectable UHE 
neutrino flux. In particular, the minimum neu- 
trino flux predicted in proton-dominated com- 
position [21] is detectable by JEM-EUSO. This 
detector can also observe the nearby UHECR 
sources using the protons with energies up to 
100 EeV and above. 

Auger clearly observes high-energy steepening 
of the spectrum, but its position and shape are 
rather different from prediction of the GZK cut- 
off. The mass composition a.t E > A EeV 
shows the dominance of nuclei becoming progres- 
sively heavier with increasing energy and prob- 
ably reaching the pure Iron composition at i? « 
35 EeV. The Auger data allows most conservative 
explanation in terms of low maximum-energy of 
acceleration, based on observation of EeV pro- 
tons. It seriously ameliorates a problem of accel- 
eration in astrophysical sources, diminishing the 
maximum energy down to (100 - 300) EeV for 
Iron nuclei. This conservative and disappoint- 
ing scenario can be conflrmed by transition from 
galactic to extragalactic cosmic rays below 1 EeV, 
by absence of cosmogenic neutrinos with the min- 
imum flux in the proton-dominated case and by 
agreement with elongation rate calculated in this 
model. 
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